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Abstract: We propose a simple higher-derivative braneworld gravity model which contains a stable 
accelerating branch, in the absence of a cosmological constant or potential, that can be used to 
describe the late time cosmic acceleration. This model has similar qualitative features to that of 
Dvali-Gabadadze-Porrati, such as the recovery of four-dimensional gravity at subhorizon scales, but 
unlike that case, the graviton zero mode is massless and there are no linearized instabilities. The 
acceleration rather is driven by bulk gravity in the form of a spin-two ghost condensate. We show that 
this model can be consistent with cosmological bounds and tests of gravity. 
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1. Introduction 

The cosmological constant problem remains one of the most challenging problems in physics. This has 
only been enhanced by the recent observational data implying the presence of a late-time dark energy 
component with w = p/p tantalizingly close to —1 Q. One line of approach to making progress with 
this problem has been to reinterpret the observed cosmic acceleration as being due to a modification of 
gravity, thus setting the true cosmological constant to zero. Although by itself this does not solve the 
problem, one may then invoke the old argument that some symmetry or selection principle may set 
A to zero. For reviews on the cosmological constant problems, see Refs. ||, [|, |), and references 
therein. 

One of the most interesting models of this kind was proposed by Dvali, Gabadadze and Porrati 
(DGP) [f[ |[ ||, [h], |ll| and has received considerable attention in recent years (see Ref. Q for a 
review). This model describes our universe as a 3-brane localized in a five dimensional bulk where 
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the brane and bulk actions are taken to be pure Einstein-Hilbert (EH) terms. As in many models 
of gravity, in which higher (or lower) scaling dimension curvature terms are included, there can exist 
multiple branches of solutions. It was shown that one branch corresponds to a self-accelerating one in 
which the brane geometry is de Sitter (dS) which could describe the observed cosmic acceleration. 

The DGP model is unusual among higher dimensional models, in that the graviton zero mode 
is massive. This implies that the graviton has five rather than the usual two propagating degrees of 
freedom and has significant implications for the physics of the model. For instance gravity can be 
shown to become nonlinear at intermediate scales |D| [l4|, |ll| |l6| . After some debate in the literature 
it has recently been demonstrated that the self-accelerating branch of the DGP model admits a ghost 
state, see Refs. jlj| H This strongly suggests that the DGP model has both a classical 
and quantum instability. 

Although a stable alternative to the DGP model has been suggested in Ref. ^0|, it relies on the 
presence of an asymmetric brane, for which the position of the brane represents a new scalar degree 
of freedom. 

In this work we propose a new model, similar in spirit to the DGP model, which also admits a 
accelerating branch in the absence of a bulk or brane cosmological constant, but for which the graviton 
zero mode is massless and no ghost states or unstable modes are present. Our model also describes a 
3-brane in a five dimensional bulk. The brane action contains an EH term and the bulk contains an 
EH and a Gauss-Bonnet (GB) term 7?. 2 . 

The key feature of our model will be to take an unconventional form for the bulk action 



where the EH term has the wrong sign. If we expand around the canonical branch vacuum, i.e. 
Minkowski space, the graviton kinetic term has the wrong sign and will lead to instabilities in the 
presence of bulk matter. By contrast if we expand around the non-canonical Anti-de Sitter (AdS) 
vacuum branch, the graviton kinetic term is conventional and no instability occurs at the linearized 
level. It is this unconventional branch that we shall consider in this work. 

This model may be understood as a type of spin-two ghost condensation pl| . Since the bulk 
EH term has the wrong sign, bulk gravitons act like ghosts, but by expanding around a non-trivial 
vacuum we recover a stable theory. In the ghost condensation work, the non-trivial solution was taken 
to be a time-dependent solution, thus breaking Lorentz invariancc p^ j. In the present case, de Sitter 
invariance on the brane is maintained by taking a solution which depends on the bulk coordinate y. 

In what follows we will: 

1. Present a summary of why the self-accelerating branch of DGP contains a ghost instability, and 
relate it with the different limits of the Randall-Sundrum (RS) model with an additional EH 
term on the brane. In particular we will show that the unstable DGP branch is related to a RS 
brane embedded in an infinite extra-dimension with warp factor diverging away from the brane, 
for which gravity is not localized. 

2. We will then propose a new model which admits an accelerating branch, which we show is stable 
at the linearized level. We emphasize that no cosmological constant is introduced on either 
the bulk or the brane. The graviton zero mode is shown to be massless and the 'dangerous' 
brane-bending modes decouple. 

3. At subhorizon scales, we show that the model effectively gives rise to four dimensional gravity 
at the linearized level. 




(1.1) 
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4. To finish, wc fine-tune the model in order to obtain the correct order of magnitude for the current 
cosmic acceleration on the brane. Ensuring that we get conventional cosmological evolution from 
nucleosynthesis onwards, requires that the five dimensional Planck mass is remarkably small 
(w 10 eV). We argue that nevertheless, this can be acceptable since gravitation is dominated 
by the brane EH term, and five-dimensional effects such as small black hole production are 
suppressed. 



2. Stability of the RS± model 

In this section we study the stability of the modified RSII model [^3| , to prepare the ground for the 
new model presented in section ^| and to gain insight into the DGP model. Readers more familiar with 
these techniques can skip straight to section |j. In the RS+ (resp. RS_) model, the extra-dimension 
has a warp factor that falls off (resp. diverge) exponentially away from the brane. In both models, we 
include a EH term on the brane: 



1 



S=—^ I d b xy/=fj R - 2A + / d^xy/^q 
2k 5 



a 1 
—^R + C m - A + —K 



(2.1) 



where R is the five-dimensional scalar curvature while R is the four-dimensional one induced on the 
brane. The bulk has a negative cosmological constant A. C m is the Lagrangian for the matter fields 
confined to the brane and K is the trace of the extrinsic curvature. The parameter a is to be associated 
to the cross-over scale r$ introduced in the DGP model, a = t§k}J 4- ^ n w hat follows, we consider 
a > 0, unless specified differently. 

When the bulk has a negative AdS curvature (RS_ model), we will show that the geometry on the 
brane is self-accelerating, but gravity is not localized on the brane, whereas when the AdS curvature 
is positive (RS+), the graviton zero mode is massless but gravity is not self-accelerating. We then 
discuss the consequences of these models when the curvature vanishes (A — > 0), and draw a parallel 
with the two branches of the DGP solution. 

2.1 Background behavior and self-accelerating solution 

In the presence of a negative cosmological constant, the bulk geometry is AdS with curvature scale \l\ 
set by A = —6/£ 2 . Since we are interested in solutions in which the brane geometry is accelerating, 
we choose the dS slicing of AdS 

d,s 2 = a 2 (y) Idy 2 + ^dx^dx") = /'f' (dy 2 - dt 2 + e 2Ht dx 2 ) , (2.2) 

sinh Hy 

where the geometry is compactified on a <S' 1 /Z2-orbifold, with a brane located at the fixed point y of 
the symmetry, y is such that sinhiJy = £H, with H the Hubble constant on the brane. The induced 
metric on the brane is hence 7^. Curved braneworlds have been considered in Refs. [^4, 25 1. 



For convenience we shall allow i to take both signs and consider only the region y > y of the Z2 
orbifold (the metric in its copy being obtained by y — -> 2y — y). If y > 0, the region is then bounded 
by the horizon at y = +00. If y < then the region is bounded by the boundary of AdS at y = 0. 
Denoting i = —e\£\, by taking c = -1 we describe the RS+ model while taking e = +1 represents the 
RS_ one. 



The extrinsic curvature is given by — — jvl + £ 2 H 2 5^ and is determined by the Israel junction 
conditions IM to be 



, „„ nl / „ l m „„\ ani / „ 1 



K$ = -fXSt -f[TZ- ^Ttt) + [R>i - -R5» ) , (2.3) 
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where T M „ is the stress-energy tensor associated with L rn . We take the canonical value for the brane 
tension: A = Q/k 2 £, and hence in our convention, the RS+ model will have a positive tension A > 
brane while RS„ has A < 0. We emphasize that due to the presence of EH term on the brane, this 
choice is arbitrary and we could have also have chosen both RS± models to have a vanishing tension 
brane, or a tension with opposite sign. This is only a matter of convention. In section |[ we shall be 
interested in a model where the brane tension vanishes, and this convention will therefore play no role. 

When some matter fields with energy density p = —Tq are confined to the brane, the junction 
conditions imply 

p/TT^ = -f<A + „ + g 



-^VTTPiP = -^(X + p) + ^H'. (2.4) 



that is, 



y/1 + PH* = 1 - e^p + ^£ 2 H 2 , (2.5) 
o 



where we write 7 = Jjjr^ > 0. We may note that in this formalism, ( 2.2) is a solution of the equations 



of motion only if the Hubble parameter is constant, i.e. only if p is a cosmological constant. However 



a more detailed calculation shows that the final Friedmann equation (2.5) remains valid for any kind 



of matter (see for instance Refs. [25, E7j for details of this derivation for arbitrary kinds of matter). 



• For the branch e = — 1, this equation admits only one solution for real H 2 . At low-energy, 
this solution is similar to the usual Friedmann equation: H 2 = K 2 p/ \£\ (3 + 27), and when no 
matter in present on the brane, H = 0. This solution is therefore the conventional one and is 
not self-accelerating. 

• For e = +1, this equation admits two branches of solutions, the canonical one for which H = 
when p = 0, and a non-canonical one for which 



when 7 < 3/2. However, if 7 satisfies this inequality, the canonical solution for which H = 
when p — couples the wrong way to matter: At low-energy H 2 = —K 2 p/ \£\ (27 — 3), and is 
not an acceptable solution. Thus there is only one physical branch for each sign of e. 

2.2 Perturbations 

To determine the stability of the model, we first consider the metric perturbations in the RS gauge 

ds 2 = a 2 (y)(dy 2 +- /ttl/ dx^dx v ) + y /Eh fiU (y,x^)dx^dx v . (2.7) 

As is usual for perturbations around a maximally symmetric spacetime, there is sufficient remaining 
gauge freedom to set /i M „ to be transverse and traceless with respect to the dS metric 7 Mt , . This is the 
statement that there are no propagating vector and scalar modes in the bulk, in the absence of bulk 
matter. In what follows we raise indices as = ^ a h av . In this gauge, the perturbed equations of 
motion in the absence of any matter in the bulk are 



;>i ■ C± ■ ^) 

y V 4 4sinh 2 J Hy 



K = 0, (2.8) 



with the notation □ = [□ — 2H 2 ] = py^V^V,/ — 2H 2 ] , where is the covariant derivative with 
respect to 7 M „. Note that in this gauge, the brane is not static but located at y = y + 5y(x^). In order 
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to derive the boundary conditions, we work instead in the Gaussian Normal (GN) frame where the 
brane is static. Performing the coordinate transformation, 



V -> V = V 



a iyY 



eff2 7 ov,u, 



the perturbed metric in GN gauge is then 



\l\ cosh Hy . 
2e\l — ,.„,„ F^dy, 



H (smhH\y |) 3 / 2 p 

where i*^,, = [V M V, y + H 2 j^, and the metric perturbation induced on the brane is therefore: 



(2.9) 
(2.10) 

(2.11) 



hnv(y) = h^ v {y) - — coth Hy F^Sy. 
In the vacuum, the perturbed boundary condition are then 



3H 

d y H coth Hy 



K(y) = -j^®K(y)- 



Translating back to RS gauge, the boundary constraint is therefore 



+ — cothiJy+^-fH 



k 2 H 



cothffyD F^y 



(2.12) 



(2.13) 



(2.14) 



Since h^ v is traceless, the position of the brane satisfies: F£ Sy = [□ + 4i/ 2 ] <5y = in the vacuum, 
or in other words, 



□ F^Sy = 2H 2 F^Sy. 



(2.15) 



Both the bulk equation (2J3) and the boundary equation ( 2.14 ) are identical for both branches e = ±1, 
but we recall, that the only difference arises in the sign of y, and in the normalization of the mode 
functions. When e = — 1, the mode functions should be normalized over the region y < y < oo, which 
requires them to fall as y — > oo or behave as plane waves, whereas for e = +1, the mode functions 
should be normalized over the region y < y < 0, which requires them to be finite as y — » _ . 

• Tensor perturbations 

We can discard for now the contribution from the brane displacement and concentrate on the tensor 
perturbations. Following the approach of Ref. |19 we write = J2m u ™(v)xffi where \ are a 
complete set of transverse traceless tensors satisfying 



2 v (m) 



Then the mode functions u m (y) satisfy 



d 2 + m 2 - 



15 



4 4sinli .Hy 

3H , „ QK? 9 

d v H coth Hy H %m 2 

J 2 y 2k\ 



H l 



u m {y) = 0, 



u m {y) = 0. 



(2.16) 

(2.17) 
(2.18) 
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It is straightforward to show that the background solution is a consistent solution of the equation of 
motion and boundary condition if m = 0: 

Mv) = A(|sinhiJy|)- 3 / 2 . (2.19) 

For e = — 1, this mode is clearly normalizable in the region y < y < oo. Consequently it describes 
the massless graviton zero mode. It can be shown that no other normalizable modes satisfy the 
boundary condition ( p. 18 ), for < m 2 < 9H 2 /A. The remaining modes are therefore massive and 



satisfy m > 9H 2 /4. They represent the continuum of Kaluza-Klein modes |2£ 



For e = +1, i.e. when y < 0, the massless mode (2.19) is no longer normalizable in the range 
> y > y since it diverges at y = 0. In fact all the modes split into normalizable modes behaving as 
(— yfl" 2 as y — > and non- normalizable modes behaving as (— y)~ 3 ^ 2 . This of course is no surprise 
since this is the familiar behavior of modes near the boundary of AdS. Demanding that the modes 
are normalizable effectively quantizes the masses m and projects out the non-normalizable massless 
mode. In appendix we show that when no matter is present on the brane, the zero mode has 
a mass m 2 > 2H 2 . Although this zero mode does not sit in the forbidden region < m 2 < 2H 2 
29 3(], |3l], |3^, |3^, |34| when the brane is empty, it will if sufficient matter is introduced on the brane. 
In other words, the brane in the RS_ model becomes unstable once the energy density confined on 
it is too important. Furthermore, since the graviton is not massless in this theory, gravity will not 
be localized on such a brane. In particular, we show in appendix [b| that at subhorizon scales, four- 
dimensional gravity is not recovered on the brane. 

• Brane-bending mode 

Let us now consider the mode which is sourced by the brane-bending term Sy. In order to understand 



this mode we follow the arguments of Refs. [|18[ [19|. The equation of motion ( 2.15 ) of the brane- 
bending mode can be associated to the one of a transverse traceless tensor with mass m 2 = 2H 2 . 
Thus although the brane-bending mode is in a sense a scalar mode, it is here degenerate with the 
tensor sector. 

Solving the bulk equation for m 2 = 2H 2 we find the general solution 

h (2) - ( A coshH V (cosh Hy — l) 2 \ pi) (9 9n x 

"* ~ { (smhHy) 3 / 2 + (sivbHyf/ 2 J ( ' 

The first solution proportional to A behaves identically to the brane-bending mode, whilst the B mode 
is distinct. 

For the RS+ model, the B mode is non-normalizable due to the divergence as y — > oo. Consequently 
any normalizable solution would have to be composed of the A mode only. This combined together 
with the brane-bending mode gives the metric perturbation in the GN gauge 




However the boundary condition (2.13) acting on this mode simply sets the combination ( Axffi — 



to zero. Thus the brane-bending decouples for the tensor sector and plays no further role. 

By contrast, the RS_ model (e = +1) requires normalizability at y — > which picks out the 
solution 

h m- Bx (2) (coshgy-i) 2 
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Then applying the boundary conditions, this mode couples together with the brane-bending mode to 
give a nontrivial solution 



U 2) 



where 



2 
H 



Z 



(coshffy- l) 2 
(sinhiiy) 3 / 2 



cosh H y 
(sinhtfy) 3 / 2 



(2.23) 



Z 



1 + ^HcothHy 



(4 + coshtfy) sinh 4 -|£ + | coth Hy + 



(2.24) 



Thus the brane-bending mode represents a genuine physical degree of freedom. This is exactly like 
the DGP scenario where the significance of this mode is explained in more detail in the literature, 
(see Refs. |l8|, |l9) for the most recent and detailed studies). When the matter content on the brane is 
such that the mass of the graviton is < m 2 < 2H 2 , the zero helicity mode of the graviton x^T°' ) is a 
ghost. In the case where m 2 > 2H 2 , xi^~°^ recovers the conventional kinetic sign (this mode is stable), 
but the brane-bending mode becomes a ghost. For the critical mass of the zero mode m 2 = 2H 2 , the 
coefficient in front of the kinetic term of x^~° vanishes and the brane-bending mode takes the role 
of the helicity-zero mode. In this case the brane-bending mode is a ghost. As long as the zero mode 
is massive, a ghost is thus present in the theory. 

2.3 Flat bulk limit and analogy with DGP 

Let us now consider the limit \£\ — > oo, (7 — ► 0), so that both the brane tension and bulk cosmological 
constant vanish which is the limit in which we get the DGP model. The stable branch of the DGP 
model corresponds to the \£\ — > 00 limit of the e = —1 RS model. It retains the same features, namely 
that the graviton zero mode is massless. The unstable self-accelerating branch corresponds to the 
\l\ — ► 00 limit of the e = +1 model. We find 



H 



2 k 2 



(2.25) 



The fact that we do not get a stable theory and a massless graviton can be seen as a consequence of the 
fact that a brane in an uncompactified AdS bulk which diverges away from the brane is pathological. 
In the limit \l\ — > 00 in which case y — > —00, the equations of motion effectively become: 



9 2 + m 2 



3H 
~2~ 



-H 2 



2^4 



u m {y) 
u m {y) 



0. 



0. 



(2.26) 
(2.27) 



This is the DGP limit and, as shown in Refs. Jll| the mass of the normalizable zero mode is 

2 2 

m = 2-^(3H — 2-^jg). When no matter is present on the brane, H = 2k 2 /o.k 2 , and the mode has 
the critical mass m 2 = 2H 2 . But as soon as some matter is present on the brane, H > 2k 2 /an 2 ,, and 
the modes lies then in the forbidden mass range < m 2 < 2H 2 which gives rise to ghosts in the scalar 
sector. 

Although the two branches of the DGP model are two natural solutions of the same model, they 
can be seen as the limit of two very distinct models |35). The stable branch of DGP e = —1, is the 
limit when \l\ — > 00 of the RS+ model, which is a well-defined stable model. The self-accelerating 
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branch e = +1 of DGP, is on the other hand the limit when \£\ — > oo of the RS_ model. Unlike in the 
RS+ case, the zero mode is not normalizable and the brane-bending mode does not decouple. This is 
the key element, responsible for the instability of the RS_ model, and by continuity, to the instability 
of the e = +1 branch of the DGP scenario. 

As a short remark, we briefly discuss here the situation when a < 0. In that case, we see that the 
self- accelerating branch is associated to the RS+ model e = — 1, where the size of the extra-dimension 

2k 2 

is finite: H = Hq = , A ■ In this case we therefore expect the normalizable zero mode to be massless, 
and hence to obtain a stable self-accelerating solution. However this model is pathological, since 
matter couples the wrong way to gravity, and this solution will therefore be unstable against matter 
on the brane. In other words there are ghosts in the brane matter sector. At low-energies, we have 

H 2 = H% &■ 



3. New Model 

In the previous sections we have seen that the presence of ghosts is intimately connected with the 
infinite volume of the extra-dimension, i.e. with the fact that the warp factor blows up away from the 
brane. In this section we propose a new model in which the extra-dimension has finite volume, and 
the warp factor falls off at infinity, but the brane geometry will be forced to accelerate in the absence 
of a cosmological constant. The key is to modify the bulk action and replace it with the action for 
a spin-two ghost condensate. To understand this, let us first consider a similar scalar field system 
defined by the action 

KW) 4 )- (3-1) 



S = J A*x^—g (+i(V0) 2 



The conventional vacuum is <fi = <po = const in which case the fluctuations 5(j) are ghosts. However, 
if we consider the scalar field to live in five-dimensional Minkowski space and perturb around the 
background solution = a y/Vd so that <j> = + 5(j), the perturbed action is then 



5,S 




Vjd^d^dcj) + (3(j 2 - l){d y 5cj>f 



(3.2) 



It is clear that provided \a\ > 1 the kinetic term will have the conventional sign and the perturbations 
will be stable. The drawback of this example is that five-dimensional Poincare invariance is broken 
down to four-dimensional. By contrast, in the gravitational version discussed below we will be able to 
perturb around an AdS background without breaking the SO{&, 2) AdS isometry group. 



3.1 Accelerating solution 

The gravitational version of the ghost condensate is defined by the action 



I 

%4 



d 5 xV=g 



-R- T *» 



d 4 x. 



a 

2k 2 



(3.3) 



where £ is an arbitrary length scale, and the dimensionless constant a is positive. The unconventional 
sign for the EH term is balanced out by the larger GB term so that perturbations around the AdS 
background solution (with length scale £) have a positive kinetic term. We stress that the GB term 
plus its boundary term only introduces terms in the action of the form {d^g vu} ) and so the equations 
of motion remain second order. In particular this means that there are only five physical fluctuations 
representing the five components of a graviton. 



- 8 - 



In section |], we will consider the limit where the five dimensional effects are highly suppressed 
k 2 3> £k\ (or in other words, when the brane is very heavy). In that limit, we need \a — 1| <C 1 in 
order to recover the right coupling of matter to gravity, but for now we leave the three free parameters 
I, k\ and a arbitrary. The GB term H2 is the trace of 

A _ t> T>A OT> A T>C orj CD t) A , r>A r> DEF 

'*~2 b ^ ^ n, c sx B — ^n, n qbd "1" def^b > 

K A BCD being the five-dimensional Riemann tensor. The boundary term Q is the generalization of the 
Gibbons-Hawking boundary terms |36|, |37j : 

Q = K + e 2 (J-2GZK»), (3.4) 

where K% is the extrinsic curvature on the brane, G^ v being the Einstein tensor on the brane and 

■\k*k%kp + \kk*kz + \ 

In the five-dimensional bulk, the Einstein equations are 



J» = --K*K$KZ + -KK^K« + -K$(K%K* - K 2 ). (3.5) 



gi = -Gi - i 2 [n 2 A B - \n 2 si j = o. (3.6) 

As already pointed out in ]3^, ^9|, at the background level, this modified Einstein equation has an 
AdS solution even though no cosmological constant is introduced in the bulk. In the dS slicing, the 
metric is expressed in (2.2). Although I could a priori be negative, we concentrate in the rest of this 
paper to the situation where I is positive: t = — e\£\, with e = — 1, using the same notation as the 
previous section. 

The boundary conditions on the brane is given by the analogue of the Israel matching conditions 
(26) in presence of GB terms @ [|| g|: 



1 „„„ \ Oi Kg ( „„ 1 

4 



n - gT*tJ -R5$ ) , (3.7) 



where 



r uvf3 — 11 otu0 



(i$ q av + R av 6% - R a p 5% - R» q a p) - -R q av - q ap 5$) , (3.8) 



i? M av p being the four-dimensional Riemann tensor induced on the brane. For cosmological implications 
on ordinary five-dimensional EH-GB models, see Refs. p2 

Working with the AdS bulk, the boundary condition ( 3.7| ) therefore reads: 

,.2 



£ 2 2 




]Vl + e 2 Hi(l + 4£ 2 H 2 ) = -^ p+ —^H 2 . (3.9) 



Provided 7 = ^j£s > 7c, (with the critical value 7c = y + 3 g 3 — 7.04), the brane geometry 
is accelerating in the absence of any stress-energy or cosmological constant on the brane. Unlike the 
previous model for which 7 was bounded 7 < 3/2, here 7 can be as large as we want. When the brane 
is empty, the geometry is indeed dS, driven by an effective cosmological constant 



3f 

M 



Aeff= 7 2T2, (3.10) 
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where r is the solution of 7 r = yjl + r (1 + 4r). For 7 > j c , this equation has two real positive 
solutions, we take f to be the smallest one. For 7> 1, this is 



A e ff 



To understand the coupling to matter, we use the notation 



11k\ pMgher = VT+PIP (1 

The Friedmann equation then takes the form 



t 2 H 2 ) - 



-£ 2 H 2 
2 



= 0(£ 4 H 4 ). 



P (7 



21 7- 



(p + Phigher) 



(3.11) 



(3.12) 



(3.13) 



where the higher order terms are negligible if £H « 1. In this limit, we recover the conventional 
coupling to matter in the Friedmann equation, if we impose the relation k 2 /2£ (7 — | ) = /tf/3, i.e. if 



3^/v^ 



1 + f +G(7- 2 )- 
2 7 



(3.14) 



When ^i? <C 1, (or equivalently, at low-energy, or when 7 3> 1), we then have 



H 2 = 



(p + A eS ) + 0(p 2 ). 



(3.15) 



We discuss the viability of this cosmological evolution in section 0. We may note again, that although 
this Friedmann equation was derived using the assumption that H was constant, (and hence for 
constant p), this result remains valid beyond this assumption and is generic for any kind of matter on 
the brane (see for instance Ref. |43| ). 

3.2 Perturbations and Stability 



Let us now consider perturbations around the bulk AdS solution in the RS gauge (2.7). This follows 
closely the discussion in section 2.2. In order to understand, how the bulk equation is affected by the 
presence of GB terms and by the overall negative sign in the bulk action, let us first consider the case 
where the bulk action is 



= d 



-9 



0£- 



^2 ( R + '—K 2 ) + £ bu 



(3.16) 



the Einstein equation gets an overall £ factor. The GB terms give an additional contribution corre- 
sponding to an extra (1 - j3) factor, (see Refs. [^jj PI or Eq.(17) of Ref.|||]). The resulting Einstein 
equation is therefore 



£(1-0) 



15 



9 



4 4 sinh 2 Hy 



5 r (U 



(3.17) 



where here again is transverse and traceless with respect to 7^ and represents the stress-energy 
of tensor matter in the bulk. Although we will consider an empty bulk in what follows, we still wrote 
down its contribution to understand the way bulk matters couple to gravity. To concentrate on our 
specific model, we consider the bulk action (3.3). This corresponds to the previous action S^'^\ where 
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now the parameters £ and (3 have been fixed to £ = — 1 and (3 = 2. The overall coefficient is therefore 
£ (1 — (3) = +1. This corresponds to a double flip of sign such that the bulk equation of motion is 
indeed unchanged and has the conventional sign: 



d 2 + B 



15 



H v 



h"„ 



(3.18) 



4 sinh iiy / 

Gravity couples the right way to matter in the bulk, even though the gravity part of the bulk action 
comes in with an unconventional sign, this is due to the presence of the GB terms that dominate over 
the bulk EH term. The bulk geometry is therefore stable against the introduction of matter in the 
bulk, at least at the linearized level, unlike the models presented in Refs. (3^, ^5| which were unstable 
due to the presence of relatively large GB terms. In this new model, the overall unconventional sign 
makes the theory stable in presence of large GB terms. In what follows, we set Tn V = 0. 

It has been pointed out in Ref. jl6| that the presence of hot matter on the brane in the case of 
an unconventional GB AdS solution leads to a naked singularity in the bulk. However our scenario 
is slightly different to that considered in Ref. [EJ, (in particular the bulk scalar curvature does not 
dominate over the GB term). It is therefore unclear whether the same result is valid here. 



Following the same procedure as section 2.2, we may work in the GN gauge (2.11) to derive the 
boundary conditions. Here again, the boundary conditions should be modified by an overall factor 



£(1 — f3) which is simply 
therefore get 



T for our specific model. Working back in terms of the RS gauge, we 



sh 

d y H — — coth Hy + 



^l-2fH cothHy □ 



M, = 



5T» - - 8T5» 
v 3 



H cothH n^t 



At H cothHy □ 



The trace of this equation imposes the equation of motion for the brane displacement: 



1 

H 



coth Hy 



U 2 H cothHy I □ I [□ + AH 2 } 5y = --^ ST. 



(3.19) 
KSy. 

(3.20) 



As mentioned in section (2.2), in the absence of tensor matter sources on the brane, Sy is such that 
[□ + 4iJ 2 ] Sy = 0. Thus the tensor F^Sy satisfies the same equation of motion ( 2.15 ) as in the 
previous model. Since y > in this model, its perturbed boundary equation is very similar to the one 
of the RS+ model which is stable. We therefore expect to recover the same features in this model. 

3.2.1 Tensors 

We focus for now on tensor perturbations in the vacuum (i.e. ST^ = 0). To start with, we do not 
consider the scalar contribution in the second line of ( 3.1 9| ) . Performing the decomposition = 

I2 m u ™{y)xffl we obtain 

'9 15 



d 2 y + m 2 - 



3H 



4 4 sinh 2 Hy 

,^2 



H' 



u m (y) = 0, 



d y + — coth Hy + — ^ - 2tH coth Hy m 



u m (y) = 0. 



(3.21) 
(3.22) 



Again it is straightforward to see that the graviton zero mode is massless m — with profile 

uo(v) = A(smhHy)-V 2 . (3.23) 

The remaining modes are massive with m 2 > 9H 2 /4. A more detailed analysis of the solutions of 
these equations can be found in appendix 
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3.2.2 Brane-bending mode 

The fact that we have found the graviton zero mode to be massless guarantees that there are no zero 
helicity graviton modes that could play the role of ghosts. Furthermore as we now explain, there is no 
brane-bending mode either and hence the modes which are dangerous in the DGP context, decouple 



in the present scenario. Following the discussion of section 2.2, in the vacuum the brane displacement 



satisfies the same equation of motion (2.15) as in the previous section. It can therefore be assimilated 



to a tensor mode with m = 2H and can thus in principle source a similar mode in the bulk. 
However, because we are interested in the same range as in the RS+ model, the normalizable mode 
with rn 2 = 2H 2 behaves identically to the brane-bending mode, i.e. u%(y) ~ coshHy/ (smhHy) 3 ^ 2 . 
Applying the boundary condition simple enforces the combination of the tensor mode and the brane 
displacement to vanish. Hence this mode does not become physical. 

Another way to argue for this is following appendix |a|, since the normalizable mode with m 2 — 2H 2 
is positive definite everywhere, it cannot be orthogonal to the zero mode which is similarly positive 
definite. But since orthogonality of different eigenmodes with different eigenvalues is guaranteed, this 
mode cannot exist. 

3.3 Response to brane matter: Modifications of Newton's law at subhorizon scales 

In what follows we study the response to a matter source ST^ on the brane. At very small scales 
compared to the Hubble radius r <C H~ , the geometry does not feel the dS expansion and appears 
Minkowski-like. Neglecting H at these scales, one can solve the mode equations exactly 



(3.24) 



where K n are the Bessel functions and the transverse and traceless part of the stress-energy is T,^ = 
8T^ V — i STjfin + 3Q<5T :M „. The induced metric perturbation on the brane h^ v in de Donder gauge is 
therefore 



h — I c T 




(3.25) 



At short-distances compared to I, i.e. for \/— Ot — > oo, the induced perturbations on the brane follow 
a behavior very similar to the usual four-dimensional one: 



/W - -^4 I ST. 



□ 



1 



:§7 



-'flis 



(3.26) 



In fact we can write this equation in the conventional form 5G M „ = k 2 ST°^, where the effective 
stress-energy is defined by 



STf v = ST^ 



(§7 



ST^u - ~ ST^ + 7^ T ;m^ 



(3.27) 



which relates in a non-local way to the stress-energy ST^. If 7 is large enough, (as will be the case 
when we put numbers in in section |J) no modification to gravity will therefore be observable at scales 
smaller that t. The fact that we recover four-dimensional gravity at small distance in presence of a 
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EH term in the brane was pointed out in Ref. p7| , (note that GB terms in the bulk have a similar 
effect). 

For subhorizon wavelengths much longer that £ i.e. £ 2 D — > 0, we have <5G M „ = n\ STff, where now 
ST% = 8Tp, + J^piy (r + log PULp, + 0(£ 4 D 2 ST^), (3.28) 

where To is the Euler's constant, (ro — 0.577). In this regime, the corrections to four-dimensional 
gravity are suppressed by £ 2 \3 and will additionally be suppressed if 7 is sufficiently large. 

The fact that we recover four-dimensional gravity on the brane at all subhorizon scales if 7 is 
sufficiently large follows from the fact that 7 measures the ratio of the brane EH term to the bulk 
EH term. When 7 is large the brane term dominates and we recover four-dimensional gravity with 
corrections suppressed by 1/7. In appendix O, we show that this remains true at the nonlinear level. 



4. Cosmological Implications 

4.1 Observed cosmic acceleration 



Our main motivation is to use the emergent acceleration presented in section 3.1 to describe the 
current cosmic acceleration. Thus the effective cosmological constant ( |3. 10j ) should be of order A e gp ~ 
f/ n\t 2 ~ 1/k\1'j I where Ijj is the current Hubble scale (size of our Universe): Ijj ~ 10 43 GeV _1 ~ 10 10 
light-years. The AdS curvature should therefore be £ ~ lnj\ff. In the regime where 7> 1, f ~ I/7. 
Our main constraint is the fact the higher orders Phigher i n the Friedmann equation ( 3.1 3| ) are 



negligible Phigher "C p i.e. p -C (7 — |) /Ik 2 , ~ 1/£ 2 k 2 , or equivalently, that £H <C 1. The initial 



theory has three free parameters «§, £ and a. As seen in section 3.1, a is fixed such that we recover 
the conventional coupling to matter in the Friedmann equation. Obtaining the correct late time 
cosmological constant will then fix the combination £k 2 . Furthermore, if we demand that the low- 
energy limit £H <C 1 or p <C 1/£ 2 k 2 is valid all the way back to nucleosynthesis (p n ~ (lO _3 GeV) ), 
then we have the additional constraint 

J?< 4 = H- 1 - l/n iy /p-^~ 10 25 GeV _1 . (4.1) 

This leads to the fine-tuning of the parameters 1/f ~ 7 ^> 1%/^ ~ 10 36 . In this limit we therefore 
have a — 1 + 9/2j ~ 1, and the constraint 7 3> 7c (necessary for the existence of an accelerating 
solution) is easily satisfied. Since 7 > 1, the five-dimensional corrections to four-dimensional gravity 
are therefore strongly suppressed. 

This model can give rise to an effective dark energy component which can be observationally 
acceptable. However we see that taking seriously the tuning of this model, the AdS length scale 
could be unusually large: I < l c = 10 3 km, i.e. of a scale at which numerous observational tests can 
be performed. One might be worried of the modification of gravity at scales of the same order of 
magnitude or below. But as pointed out in ( |3.27 ), if 7 is large enough (and in our model it could be 



at least 10 36 ), the modifications of gravity at subhorizon scales will be unobservable. 
4.2 Five-dimensional Planck scale 

The low-energy constraints leads to a five-dimensional gravitational coupling constant k 2 3> k 2 , where 
k 2 = ljjK 2 /a£ c ~ (lO _8 GeV) . This corresponds to a remarkably low five-dimensional Planck mass 
of about 10 eV. The fact that the five-dimensional Planck scale is so small, suggests that we would be 
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able to easily produce small black holes in current laboratory experiments, something which is clearly 
ruled out. However we argue that this is not necessarily the case. Consider the collision of two high 
energy particles on the brane, the total stress-energy of the brane is not given just by the stress-energy 
of matter, but by the combination 



4 



As is clear from the arguments given so far, (e.g. equation ( 3.27 )) these two contributions cancel out 
at leading order and give rise to 

Tup ^ T, L i/ ^ n~Tnv. (4-3) 

7 4 

Thus for large 7 the effective stress-energy on the brane is a negligible fraction of the real matter 
stress-energy. Indeed it is suppressed by the same factor of 1/k 2 , as the bulk action. Thus the system 
behaves like the RS model but with an effective five-dimensional gravitational coupling k 2 , — k 2 £ which 
corresponds to a five-dimensional Planck mass of an acceptable M5 = 10 4 — 10 5 GeV. We expect then 
that the production rate of black holes will be determined by this scale, which is sufficiently high 
to have evaded current accelerator experiments. In appendix C we give some nonlinear solutions 
describing high energy particles which exhibit this effect. If additional dimensions are compactified 
on very small scales the fundamental higher-dimensional Planck scale could again be much larger. 

4.3 Finite maximal energy density 

To finish, we remark briefly, on the "maximal energy density" feature generic to RS-GB models, as 



pointed out in Refs. |4<j| . Considering the Friedmann equation (3J3), there exist a maximal finite 
energy density 



7 3 , 7 



O(^), (4.4) 



for which the Hubble parameter is H 2 ~ 7 2 /^ 2 . If this model was trusted all the way back to the Big 
Bang time, the Big Bang temperature, scale factor and Hubble constant would be finite. The solution 
is nevertheless singular since H diverges, this is a very mild unusual spacelike singularity. Although this 
maximal value would be huge, (p max ~ (10 15 GeV) 4 with the previous fine-tuning), it is nevertheless 
comparable to the energy scale typically taken at the beginning of inflation (pbcg ~ (10 16 GeV) 4 , see 
Ref. which suggests that the inflation should start at much lower energy in this model. The 

end of inflation is usually taken to be at around 10 13 GeV, so this model could strongly affect the 
behavior of some inflationary models. The study of this feature is however beyond the scope of this 
paper and we suggest Ref. E8| for further details on this interesting issue. 



5. Perspectives 

We have presented a model of gravity which accelerates in the absence of any cosmological constant or 
potential which can provide an equivalent description of most of cosmic history without the need for 
the introduction of a cosmological constant. There are many questions that this model raises. Firstly, 
the formal: Is the bulk theory stable beyond linear order? If not what do the instabilities imply? 
Are there other higher derivative actions with similar properties? In particular, can we extend it to 
models where the bulk action takes the conventional sign but with higher order corrections? It has 
recently been demonstrated that higher derivative corrections to gravity alone, to all orders in the a' 
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expansion, give rise to AdS solutions in string theory, even in the absence of a cosmological constant 
term induced from the matter sector, see Ref. [^o|. Our model could be viewed as a toy version of 
these. Can such an action be derived from a UV complete theory? Secondly, the phenomenological: 
How do the high energy modifications of gravity affect inflation, structure growth? How well can we 
constrain the free parameter £? And finally how can the different cosmological tests performed on the 
DGP model (see for instance Refs. 51 1) be extended to this scenario? 
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A. Eigenvalues and eigenfunctions for the new model 

Let us now analyse the normalizable tensor perturbations for the new model. To simplify the equations 
of motion and boundary conditions we define /j, = m/H, x = Hy. The brane position is therefore at 
sinh 5 = £H. The the bulk equation then becomes 



dz 



9 



15 



4 sinh 



u m (x) = 0, 



(A.l) 



and the boundary condition is 



d x 



coth x 



■ coth x(l — 2 sinh 2 x] 



u m {x) = 0. 



(A.2) 



Here we have made explicit use of the background equations of motion when the brane is empty, i. e. 
Eq. (|3.9[). The inner product is defined as 



(/,<?) 



I 

dxf(x)g(x) + — coth 5(1 — 2 sinh x)f{x)g{x) 



(A.3) 



and one can demonstrate that non-degenerate eigenfunctions will be orthogonal with respect to this 
definition. To ensure that the norm is positive definite we require that sinh a; < \j\[2 which is the 
case in our model. For /i 2 > 9/4 there are a continuum of normalizable modes which behave like a 
superposition of plane waves at the horizon. In the range < fi 2 < 9/4 normalizability requires that 
the modes fall of as exp(— fix) as x — > oo where p, = \/9/4 — /i 2 . The solution Eq. (A.l) with this 
property is 



t (x) — A M v sinh a; (cosh a;) , 



(A.4) 



where Cffl is the associated Legendre function and v — — | + A* 2 - The existence of a mode in the 

region depends on whether this function satisfies the boundary condition (A.2). It is straightforward 
to show that for /i 2 = the solution is 



u (x) 



A 
(sinh x) 3 / 2 



(A.5) 
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and that this is normalizable and consistent with the boundary condition. Now for any normalizable 
function for /i 2 < 9/4 we have d 2 u/u > and as x — > oo, lnu — > — \/9/4 — m 2 . Assuming u > 
the gradient c^w increases with increasing x, which in turn implies c^it < for all x and hence the 
functions A. 4 have no zeros, i.e. there is no real solution to u(x) = 0. Since the definition of the inner 
product is positive definite, it would be impossible for such a mode to be orthogonal to the zero mode, 
and hence the the only normalizable mode consistent with the boundary conditions in the entire range 
/j, 2 < 9/4 is the zero mode. This result is crucial because it implies both that there are no ghost states 
arising from the dangerous regime < fi 2 < 2 and no unstable modes with /j 2 < 0. 



B. Gravity in the RS model 



In this section, we evaluate the mass of the zero mode in the RS_ model of section 2.2. Using the 



same notation as the previous appendix, the bulk equation of motion is the same as (A.l), and for 
jj 2 < 9/4, the solution of the bulk equation which is normalizable in the range x < x < is 



,(x) = A^\J sinhxPj 2 ) (coshx) , 



(B.l) 



where is the other associated Legendre function. This solution is valid if /i 2 ^ and /i 2 ^ 2. For 
/.t 2 = 0, the normalizable mode in the region x < x < is 



u (x) = A Q 
and for /i = 2, the normalizable mode is 

«2 = A : 



cosh x — 3 cosh x + 2 



(sinh x 



3/2 



(coshx — 1) 
(sinhx) 3 ^ 2 



(B.2) 



(B.3) 



If we consider the brane to be empty, p = 0, the Hubble parameter is given in (2.6) and the boundary 
condition is therefore 



3 / 2 

d x + - coth i + /iVl--7 



u m (x) = 0, 



(B.4) 



where sinhx = ~^yl — §7- We can easily check that both (B.2) and ( [B.3| ) never satisfy this boundary 
condition no matter what 7 is. However, when ji 2 = 2, the mode can couple to the brane-bending 
mode. Effectively , the presence of a brane-bending mode with /x 2 = 2, appears as a source term in 



the boundary condition (B.4). This possibility is explored in section 2J2 and is shown to be the source 
of an instability on this brane. 



Using the more general solution (B.l) on the other hand, we see that the boundary condition (B.4) 
is satisfied for /1 2 = n 2 (j) > 2, where its exact value depends on 7. In this model, despite the fact 
that the zero mode is massive, it does not lie in the forbidden region < /1 2 < 2, provided no matter 



is introduced on the brane. For H greater than the vacuum value, (2.6), the value of pi 2 (j) decreases 
and crosses 2 when H is greater than a critical value i.e. when enough matter is introduced on the 
brane. The model becomes then unstable. Although this self-accelerating branch appears stable under 
a critical value of the energy density on the brane, we can show in what follows that gravity cannot 
be localized on the brane. 
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At subhorizon scales, the Hubble parameter is negligible, and the equations of motion become 



„ 3 ani 



K{y) 



^b 1 vi 



(B.5) 
(B.6) 



where y — £ < and T^ v is the stress-energy tensor introduced on the brane, which is assumed to be 
traceless for simplicity. The normalizable mode function in the range £ < y < is h ~ \fyl 2 {^/— Dy) 
where /„ is the Bessel I function. On the brane, the mode functions are therefore: 



KAv) 



-□) 



Ik 



■-Uh 



-T 



so that the perturbed metric on the brane is 







\ +3 2^jf n + ... 


_ 4 







-□) 



(B.7) 



(B.8) 



Thus the induced perturbations on the brane don't satisfy the ordinary four-dimensional equation, 
and gravity will not appear localized on this negative tension brane. This can be understood easily 
by the fact that the warping factor increases exponentially away from the brane and gravity will tend 
to be localized at the boundary of AdS y = 0, i.e. infinitively far away from the brane, unlike the 
positive tension brane case. 

2 

Writing ~ ro, the cross over-scale, one however recovers the standard DGP result in the limit 

\£\ -> oo 



-□ - r □ 



(B.9) 



such that at very small scales k » l/Vo, one recovers the four-dimensional standard behavior h^ u 

2 

— -^T^v, whereas at large cosmological scales, gravity becomes five-dimensional. 



C. Nonlinear gravitational waves 

An ultra-relativistic particle may be described by a brane stress-energy of the form T uu = E8 2 (x)5(u) 1 
where E is its energy. For £ -C H^ 1 we can find an exact nonlinear solution of the bulk equations in 
the form of a warped pp-wave 

ds 2 = a 2 (y) (dy 2 + 2dudv + a _3/2 (y)/(y, u, x 2 , x 3 )du 2 + dxj + dxfj , (C.l) 

where u = (xi — t)/V2, v = (xi + t)/V2 and a(y) = l/y. This is a solution provided both the bulk 
equation and the boundary condition 

(C.2) 

k, 2 T uu (C.3) 



f(y,u,x 2 ,x 3 ) = 



2£ □ 



f(y,u,x 2 ,x 3 ) = - 
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are satisfied. Here □ = d 2 2 



dl 3 



and the brane is located at y = y — £. In fact this is identical to the 
special case of linearized perturbations considered earlier, except the solution is now valid non-linearly 
for / of arbitrarily large amplitude. The general solution for arbitrary source T uu is 



/ = - 



34 



£( 7 -3) 




K (h 



2) t^/^UK {h 



It 

□ " 



Similarly as in section 3.3, at very small scales, £ 2 D 3> 1, this simplifies to: 

3(7-1) 



/(r<^) 



(7-|) (7-6) 

whereas at scales much larger than £, / simplifies to 



b_rp 

— , 11. 



f(r>i) ~ - 



3«§ 



£( 7 -3) 



1 

□ 



2(f -1) 



r + log : 



(C.4) 



(C.5) 



(C.6) 



The solution corresponding to T uu = E 5 2 (x)5(u) is the five-dimensional version of the Aichelberg-Sexl 
metric |32| generalized to the present model. Similar shock-wave solutions in the DGP context have 
been considered in |53, M. In particular at small scales, / is 



3(7 



(7- 



(7 



9\ 2 
2 ) 



(7- 



!)(7-6)^ )l0 4 



(C.7) 



where r = \/a 



At small scales, we therefore recover the standard four-dimensional result, up 



to a factor which tends to 1 when 7 ^> 1. 

At larger scales, when H^ 1 ^> r ^> I, the corrections to the standard result are suppressed by a 



factor of y- < 1: 



f(r>i) 



k\ (7 



1" (7- 3) 



log- 



f 2 



2(f 



l) r 



-A 



ES(u), 



(C.8) 



where A is the cut-off scale. Typically in this approximation, we can take it to be A ~ r/£, so that 
the correction term is suppressed by both a factor of 1/7 and a factor of £/r. 

As long as 7 > 1 this solution is essentially equivalent to the conventional four-dimensional 
Aichelberg-Sexl metric on the brane. Since this stress-energy source is traceless, the position of the 
brane remains unperturbed by the gravitational wave. This solution differs significantly from that say 
in the RS model, where for distances much less than £ we recover the five-dimensional version of the 
Aichelberg-Sexl metric. This confirms our general argument that the production rate of black holes 
from high energy collisions is suppressed. 
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